A Search for the Electric Dipole Moment of the r-Lepton 



The ARGUS Collaboration 

H. Albrccht, T. Hamachcr, R. P. Hofmann, T. Kirchhoff, R. Mankel\ A. Nau, S. Nowak^, 
D. Refiing, H. Schroder, H. D. Schulz, M. Walter^, R. Wurth 

DESY, Hamburg, Germany 

C. Hast, H. Kapitza, H. Kolanoski\ A. Kosche, A. Lange, A. Lindner, M. Schieber, 
T. Siegmund, H. Thurn, D. Topfer, D. Wegener 
Institut fiir Physik^, Universitdt Dortmund, Germany 

C. Frankl, J. Graf, M. Schmidtler* M. Schramm, K. R. Schubert, 
R. Schwierz, B. Spaan, R. Waldi 
Institut fiir Kern- und Teilchenphysik^, Technische Universitdt Dresden, Germany 

K. Reim, H. Wegener 
Physikalisches Institut^, Universitdt Erlangen-Nitrnberg, Germany 

R. Eckmann, H. Kuipers, O. Mai, R. Mundt, T. Oest, R. Reiner, A. Rohde, 

W. Schmidt-Parzefall 
//. Institut fiir Experimentalphysik, Universitdt Hamburg, Germany 

J. Stiewe, S. Werner 
Institut fiir Hochenergiephysik'^ , Universitdt Heidelberg, Germany 

K. Ehret, W. Hofmann, A. Hiipper, K. T. Knopfle, J. Spengler 
Max-Planck-Institut fiir Kernphysik, Heidelberg, Germany 

P. Krieger*, D. B. MacFarlane^, J. D. Prentice®, P. R. B. SauU^, K. Tzamariudaki^, 
R. G. Van de Water®, T.-S. Yoon® 
Institute of Particle Physics^'^, Ganada 

M. Schneider, S. Weseler 
Institut fiir Experimentelle Kernphysik^^, Universitdt Karlsruhe, Germany 

M. Bracko, G. Kernel, P. Krizan, E. Kriznic, G. Medin^^, T. Podobnik, T. Zivko 
Institut J. Stefan and Oddelek za fiziko^^, Univerza v Ljubljani, Ljubljana, Slovenia 

V. Balagura, S. Barsuk, I. Belyaev, R. Chistov, M. Danilov, V. Eiges, L. Gershtein, 
Yu. Gershtein, A. Golutvin, O. Igonkina, I. Korolko, G. Kostina, D. Litvintsev, P. Pakhlov, 
S. Semcnov, A. Snizhko, I. Tichomirov, Yu. Zaitsev 
Institute of Theoretical and Experimental Physics^^, Moscow, Russia 



^ Now at Institut fiir Physik, Humboldt-Universitat zu Berlin, Germany. 



^ DESY, IfH Zeuthen, Germany. 
Supported by the German Bundesministerium fiir Forschung und Technologie, under contract number 054DO51P. 
Now at Caltech, Pasadena, USA. 

Supported by the German Bundesministerium fiir Forschung und Technologie, under contract number 056DD11P. 
^ Supported by the German Bundesministerium fiir Forschung und Technologie, under contract number 054ER12P. 

Supported by the German Bundesministerium fiir Forschung und Technologie, under contract number 055HD21P. 
® University of Toronto, Toronto, Ontario, Canada. 
^ McGill University, Montreal, Quebec, Canada. 

Supported by the Natural Sciences and Engineering Research Council, Canada. 

Supported by the German Bundesministerium fiir Forschung und Technologie, under contract number 055KA11P. 
On leave from University of Montenegro, Yugoslavia 
^"^ Supported by the Ministry of Science and Technology of the Republic of Slovenia and the Internationales Biiro KfA, 
Jiilich. 

^'^ Partially supported by Grant MSB300 from the International Science Foundation. 



Abstract 



Using the ARGUS detector at the e^e~ storage ring DORIS II, we have searched for the real and 
imaginary part of the electric dipole formfactor dr of the r lepton in the production of r pairs at 
= 100 GeV^. This is the first direct measurement of this CV violating formfactor. We applied the 
method of optimised observables which takes into account all available information on the observed 
r decay products. No evidence for CV violation was found, and we derive the following results: 
7^e(d^) = (1.6 ± 1.9) • 10"^^ecm and Im{d^) = (-0.2 ± 0.8) ■ lO^^^ecm, where statistical and 
systematic errors have been combined. 



The understanding of CV violation and the question if the current description of the elcctrovifeak inter- 
action in the mimimal Standard Model (MSM) is sufficient to explain the matter-antimatter asymmetry 
in the universe is one of the burning issues of present particle physics. It appears that there is insuffi- 
cient CV violation in the MSM to generate the observed baryon asymmetry [Q, stimulating the search 
for additional sources of CV violation beyond the MSM, e. g. leptoquarks, additional Higgs bosons or 
supersymmetric particles. 

One way to search for new CV violating interactions is in the pair production of r leptons. Contributions 
from such interactions are parametrised model independently by the electric dipole formfactor dr{q^). 
Up to now, upper limits on the weak dipole formfactor have been determined at LEP ||^, and indirect 
measurements exist of the electric dipole formfactor Q . 

The T dipole formfactor is measured in charge dependent momentum correlations in the final states of 
the reaction e^e~ t~^t~ f^f~- At low energies, r pairs are produced in the CV even state ^5*1 via 
a virtual photon with quantum numbers J^*^ = 1^^. CV violation would result in the CV odd state ^Pi 
which leads to other correlations between p{f~^) and p{f~) than in the state '^S'l. 

In the Born approximation, the Lorentz invariant matrix element for the production of r pairs in the 



e — > 7 



T+T including CP-violating terms is given by ^ 



lA^prodl OC ipiq2){P2qi) + {Piqi){P2q2) + m^{piP2) 

+ (wi^a) (^{piq2){P2q2) - ^ml{piP2) 

- {uJlPl){uJlP2){P2q2) - (wiPl)(w2P2)(Piq'2) 
-I- ne{dr) (2 rUr e^^^pa [Pl +P2T^^{i^2 

{Pi {{P2qi) - (piqi)) - qi{piP2))) 

+ Im{dr) {{ml{uji - uj2)f, {-{piqi)Pi + {P2qi)Pi + (pi9i)P2 " (P2'7i)P2) 
+ ('Zi'72) ((pigi)(-(t^ipi) + {uj2Pi) + {(^iP2) - (W2P2)) 

+ {P2qi)iiuJlPl) - {i02Pl) - {l^lP2) + [t^2P2))) 

- {0J2qi) {{piqif + {P2qif) + {uJiq2) {{piqif + {P2qif) 
+ 2 {piqi){p2qi) {{(^2qi) + {uJiq2)))) 
+ {-ml{pip2) - {pm)^ + 2{piqi){p2qi) - {p2qif + (piP2)(gi'72)) 

{ml + {qiq2) - (t^2<7i)(t^ig2) + ml(uJiL02) + (gi'72)(t^iW2)) , 

where (iT-(s) is the electric dipole formfactor of the r lepton, pi_2 {qi,2) are the 4-momenta of e^,e~ 
(t+,t^), and uji^2 denote the 4-polarisations of t+ and r^. 

Because of the long lifetime of the r-lepton, the matrix elements for the production of a r pair, |AlprodP, 
and that of the subsequent decays, |A^docP, factorize in the Born approximation: 

= |Xp,odnXdcc|'. 

l-A^prodP shows a linear dependance on TZe{dr), Tm{dr), and Id^-p, whereas |AldccP is independent of the 
electric dipole moment. Therefore, the differential cross section of the process e+e~ — > r+r^ f~^f~ 
will exhibit a similar linear dependance : 



datot = dasM + 'Re{dr)dan + Jm{dr)dax + \dr\'^daQ, 



where dusM corresponds to the differential cross section as expected in the Standard Model. Since 
\dr\'^ • s ^ 1, IdTpdcTQ can be neglected for \dr\ < 5 • 10~^^ ecm, and the differential cross section can be 
rewritten as: 



datot = dasM ■ [1 + Tle{dr)uin + 'lT[v{dr)uji\ , lo-ji = , ujx = ^ 



dasM ' dasM 

In this form, uj-ji and ajj are the so-called optimised observables as suggested in ref. [Q . Both quantities uj-ji 
and uji describe CV asymmetries with the properties / ajjidasM — J Aiders m — and J uj-jzuJxdasM — 0. 
This leads to 

T^{dr) = — , Jm(d^) = —2- , 
where the means are defined by 

r LOj datot 
^^^^ = -Jd^ ' 

with j = TZ,T, and (w^) accordingly. For the statistical errors, error propagation gives: 



2\ 



Experimentally, the only accessible means are 

_ ^ujj 1] datot 
J datot ' 

(w|)cxp accordingly, and real and imaginary part of c?r are determined by (a;j)oxp/(w|)oxp- The influence 
of the detector acceptance rj has to be determined by Monte Carlo simulation and will be discussed later 
in this paper. 

The following analysis is limited to the final states tt [pi'){pv)^ {pi'){ijLVv) and {pi'){evv) since these 
are the most sensitive ones given the ARGUS environment. The application of the method of optimised 
observables requires the integration of each observed event over all unmeasured quantities, such as r 
direction, photons of the initial state bremsstrahlung which mostly escape undetected, radiated photons 
in the decay l^vv that merge in the cluster of the charged particle and photons from external 

bremsstrahlung. 

The integration is performed numerically using Monte Carlo methods. 500 tries are made for each event 
using a hit and miss approach for the relevant kinematical quantities. As in refs. Q and Q], where this 
method has been used for the determination of Michel parameters, there are rihit successful tries in which 
the observables are compatible with the kinematics of a rr ^ 71/1/2(72) event. The resulting differential 
cross section for each observed event is then the mean of these nhit values: 



da(dr) = daAdr) 



"hit 

The validity of the method for determining dr has been tested by generating events with the KO- 

RALB/TAUOLA Monte Carlo program [||. To describe the effect of an electric dipole formfactor, these 

events were weighted with the full matrix element Q . Details can be found in ref. ^ . 

The analysed data sample has been collected between 1983 and 1989 with the ARGUS detector at the 

e"'"e~-storage ring DORIS II. The detector and its trigger requirements are described in ref. The 

integrated luminosity for this analysis was 291 pb~^ corresponding to 277500 r pairs produced. 

The pp event selection follows closely the one presented in ref. |lO[ and starts with requiring two oppositely 

charged tracks forming a vertex in the interaction region and with —0.9 < cos(piP2) < —0.5. This reduces 



the background from more collinear QED events, like Bhabha and /z pair events, and takes into account 
that the deacy products of the two t are in opposite hemispheres. Each charged track must have a 
transverse momentum above 60MeV/c and point into the barrel region (| cos9\ < 0.7) to ensure good 
trigger conditions. Neutral pions are reconstructed by their 77 decays. Photons are required to deposit 
at least 100 MeV in the calorimeter. The number of photons in the event has to be 2, 3, or 4. Each of 
these photons must have an angle between 10° and 90° to one of the two tracks. The 10° cut is against 
faked photons from charged particle split-offs. The photons are then assigned to their closest track. Only 
one or two photons are allowed per track. Depending on the number of detected photons, neutral pions 
were reconstructed in two ways. If only one photon was found, a minimum energy of 1 GeV was required 
for this photon to be taken as single-cluster-iT^ . In case of two photons the system had to fulfill < 9 
when kinematically constrained to the 7r° mass. An additional cut was made on the combination of the 
charged track and the reconstructed 7r° to fulfill m^±^o < 1.8GeV/c^. Two-photon and QED reactions 
are suppressed by the following requirement: 

l5IPT|c/Scms > (0.45 •^|plc/£;e„,s- 0.5)2 + 0.01) ^ 

where p denotes the momentum and px denotes the transverse momentum of charged and neutral parti- 
cles. 

The pi selection was performed like in ref. 0. Wc only present here the differences to the pp selection. 
The opening angle of the charged tracks must fulfill —0.965 < cos(pip2) < 0.2. This cut is much weaker 
because the background of hadronic events is much less in this decay channel. We also use two different 
I ^ Pt I requirements depending on the tt*' reconstruction: 

single cluster : l^pVlc/Scms > 5 ■ {M^^^J E^^, - Ab^ , 
two cluster: |^PT|c/i?cms > 3 • (M^ j^JiJ^^, - 0.4)^ , 

where M^j^g ~ (-Ecms ~ -^i)^ ^ (X^Pi)^- To ensure a good lepton identification, a momentum above 
0.5GeV/c for electrons and above l.lGeV/c for muons was required. A possible lepton candidate had 
to fulfill the following particle identification cuts: for muons > 0.8, Lg < 0.2, and for electrons 
Le > 0.8, Lfj^ < 0.2, where the likelihood ratios Le.^i are defined in ref. In Fig. |l|, the resulting tttt mass 
distributions for accepted events are compared to Monte Carlo calculations for pure signal events. The 
differences between those spectra are mainly due to r background from the decays tt (pi^)(7r7r°7r*'i^) 
and TT —> {lh'i'){7nr^TT'^v). 
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Table 1: Monte Carlo results for the background estimation. 

The distribution of nhit is shown in Fig. ||. The good agreement between iV(nhit) for selected data and 
accepted Monte Carlo events supports the validity of the assumptions used. The high value of (rihit) 
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Figure 1: 

TT TT mass spectra for selected (a) (p i>){p^i/) and (b) (/x events. The points 

with error bars represent the data. The hatched histograms show the expectation from the 
KORALB/TAUOLA Monte Carlo. 



shows the effectiveness of the numerical integration. A cut of rihit > 35 is applied. After the selection we 
get 1074 {pv){pv), 2052 {pi'){ij,i'i/), and 3675 {pv){evv) candidates at a mean center of mass energy of 
-y/i = 10.38 GeV, defining the scale for the measurement of dr. 
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Figure 2: 

Distribution of rihit for (a) [p v){p^i') and (b) [p vh'){p^i') where rihit is the number of 
possible kinematic configurations in 500 tries. The points with error bars represent the 
selected data. The hatched histograms show the expectations from KORALB/TAUOLA. A 
cut of rihit > 35 is applied, lower rihit values originate mostly from events with hard initial 
state radiation. 

The selected events lead to distributions of the optimised observables as shown in Figs. ^, ^, and ||. 
The first and second moment of these distributions are used to determine the electric dipolc formfactor. 
All six distributions are centered around zero indicating dr close to zero. The widths of the three 
distributions are very similar. Note that the width of the uj-n distribution is largest for the pp channel. 
This reflects the higher information on Tie{dr) per event. 

The results have still to be corrected for systematic uncertainties, like acceptance corrections, particle 
identification, and background. The influence of background channels was studied using Monte Carlo 
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Figure 3: Optimised observables in the decay rr {pv){piy) 
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Figure 4: Optimised observables in tlie decay rr — > {pv){yLw) 



calculations. For the {pv)[evv) and {pv)[^vu) channels the background of hadronic, gamma- gamma and 
QED events (Bhabha, /x-pair) was found to be negligible. The only background source is from non-signal 
T events. In the {pv){pv) channel, the background from gamma-gamma and hadronic events has to be 
considered. The resulting background contributions are summarised in Table |l|. 

The influence of the r background on dr was studied with event samples generated with dr — 0. This is 
possible because of the small amount of background and the loss of spin information in most channels, 
e. g. in the decay — > ■k^-k^tt^v where some photons or a tt" were not detected. The Monte Carlo 
results were used to correct the obtained values for 7?e((i^) and Im{dr). Results after correction are 
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Figure 5: Optimised observables in the decay rr {pi/){ei'u) 



given in Table 3 with their statistical errors. 



systematic error source 


real part [10 ^®ecm] 

KPP J)Pt^ J)P'' 
"TZc "TZc "TZc 


imaginary part [10 ^^ecm] 
j)PP /;PP AP" 


acceptance correction 
charged asymmetrie 
lepton identifikation 
photon identification 
background variation 


0.48 1.00 1.00 

0.022 0.010 
0.005 0.002 0.002 
0.06 0.13 0.001 


0.35 0.38 0.36 
0.24 0.06 
0.007 0.006 
0.003 0.002 0.001 
0.02 0.02 0.003 


total systematic error 


0.48 1.01 1.00 


0.35 0.45 0.37 



Table 2: Contributions to the systematic error. 



Detector acceptance effects have been studied in a Monte Carlo simulation by generating large event 
samples with Tie{dr) and Im{dr) equal to 1 • 10^^® and 5 • 10~^^ecm. After reconstruction including 
all acceptance cuts, these values were obtained again within uncertainties listed as systematic error in 
the first line of Table 3. Strong infiuences on the dr measurement may occur if the detector acceptance 
is different for particles of opposite charge. Therefore, this effect was extensively studied with simulated 
data where holes in momentum, polar and azimuthal angle were generated. No charge dependent effects 
were found. In addition, systematic errors in the determination of the momentum of charged particles 
showed no influence on the dr determination. A detailed study was necessary because charge dependent 
systematic momentum shifts of up to ±20 MeV/c for momenta above 4 GeV/c have been observed for 
some data taking periods. 

The systematic error estimates are given in Table |[ their sums also in Table ^. No evidence for CV 
violation was found. One of the three final states {pvpv) was also investigated with a maximum likelihood 
fit method. The result of this fit was identical to the one with optimised observables, the statistical errors 
were equal within 1%. 



channel 


dipolc formfactor 


value [10 ^^ecm] 


PP 


TZe{dr) = 3.3 ±2.5 ±0.5 
Jm(d^) ^ 0.5 ±1.4 ±0.4 


PP 


7^e(^^^) = -4.0 ±5.3 ±1.0 
Xm(dr) = -1.2 ±1.8 ±0.5 


pe 


7^e(d^) = 0.2 ±3.1 ±1.0 
Xm(dr) = -0.2 ±1.1 ±0.4 



Table 3: Final results with statistical and systematic errors for all three selected channels. 



With statistical and systematic errors added in quadrature, the combination of the three measurements 
leads to: 

7^e(^^^) ( 
Xm(dr) = 

which translates into the following upper limits 

|7^e(d^)| < 4.6- 10"^^ ecm (95%CL) , 
|Xm(d^)| < 1.8-10"i^ecm (95%CL) . 

The result is in agreement with current LEP results which used the assumption of a non-existing 
anomalous magnetic dipole moment to calculate an upper limit on the electric dipole moment. 
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